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ABSTRACT 


This  paper  describes  aircraft  observations  made  with  the  NOAA  P-3  research  air¬ 
craft  on  24  March  1989.  The  measurement  region  was  over  the  marginal  ice  zone, 
southeast  of  Spitzbergen,  between  74  to  76°  N  and  19  to  27°  E.  Above  the  atmospheric 
boundary  layer,  the  geostrophic  wind  was  12  ms-1  at  approximately  15°  ofT-ice.  A  well 
mixed  layer  extended  along  the  entire  200  km  flight  path  which  was  perpendicular  to  and 
centered  over  the  ice  edge.  There  was  a  stratocumulus  layer  over  the  ocean  which  de¬ 
creased  in  thickness  towards  the  ice.  The  associated  inversion,  which  coincided  with  the 
cloud  layer  top,  continued  to  decrease  in  height  over  the  ice.  This  caused  a  strong 
thermal  wind  effect  within  the  atmospheric  boundary  layer.  The  near-surface 
geostrophic  wind  decreased  from  approximately  12  ms-1  over  the  ocean  to  5  ms~'  over 
the  ice  due  to  the  horizontal  temperature  gradient  and  sloping  inversion. 
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I.  INTRODUCTION 


A.  ARCTIC  ICE  EDGE  METEOROLOGY 

Desire  to  make  use  of  the  Arctic  polar  region  resources  has  greatly  increased  man's 
acti\ity  there.  Understanding  the  environment  is  vital  in  order  to  tap  the  immense  dis¬ 
coveries  of  oil,  natural  gas  and  rich  fisheries.  The  relationship  between  global  weather 
and  the  polar  ice  cap  has  also  been  an  item  of  long  standing  interest.  In  the  Arctic, 
between  1953  and  1977,  the  average  maximum  extent  of  the  sea  ice  was  14  million  km2, 
which  is  5%  of  the  northern  hemisphere.  The  minimum  extent  was  7  million  km2  (Walsh 
and  Johnson,  1979).  How  this  annual  variation  in  ice  cover  is  related  to  large-scale  and 
mesoscale  circulations  is  not  fully  understood. 

The,  marginal  ice  zone  (MIZ)  is  the  transition  region  between  the  open  ocean  where 
there  is  no  sc  a  ice  and  the  polar  ice  pack  with  80%  or  greater  ice  cover.  Large  horizontal 
gradients  in  ice,  ocean  and  atmospheric  properties  characterize  the  region.  Also,  strong 
interactions  between  the  ice,  ocean  and  atmosphere  occur.  The  spatial  variation  of  the 
atmospheric  boundary  layer  (ABL)  thermal  structure  and  surface  roughness  differences 
between  the  ice  and  ocean  cause  secondary  atmospheric  and  oceanic  circulations  near 
the  icc  edge. 

Warm  ocean  currents  limit  the  extent  of  the  ice  cover.  The  wind  controls  the  con¬ 
centration  of  icc  in  the  MIZ  (Davidson  and  Guest,  1987).  The  ice  characteristics  from 
the  polar  ice  cap  across  the  MIZ  vary.  Overland  et  al.  (19S3)  describe  the  MIZ  by  three 
ice  regions:  1)  a  smooth,  interior  ice  pack  region  of  large  ice  floes  with  1  km  horizontal 
and  2  to  3  m  vertical  dimension;  2)  a  transition  zone  of  small  rectangular  floes  with  20 
m  horizontal  anJ  3  to  6  m  \ertical  dimension;  and  3)  a  rough  zone  of  small,  broken  floes 
with  less  ice  concentration  than  zone  two.  These  varying  conditions  of  surface 
roughness  affect  the  wind  and  interactions  between  the  atmosphere,  ice  and  ocean. 

The  ice  is  relatively  fresher  than  the  ocean.  The  melting  and  freezing  of  the  ice  can 
affect  the  temperature,  salinity,  and  vertical  density  structure  of  the  upper  ocean.  The 
change  in  surface  stress  across  the  MIZ  is  belie\cd  to  cause  ocean  upwclling  near  the  ice 
edge  (McPhee,  I9S3). 

The  atmosphere's  temperature  markedly  varies  across  relatively  short  distances, 
from  warmer  temperatures  over  the  ocean  to  very  cold  over  the  ice.  Atmospheric  sta¬ 
bility  can  also  experience  abrupt  variation.  Over  the  ocean,  the  atmosphere  is  typically 
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unstable,  while  over  the  ice  it  is  stable  to  slightly  unstable.  Strong  inversions  are  com¬ 
mon.  Shapiro  and  Fedor  (1989)  observed  very  intense  cyclones,  referred  to  as  polar 
lows,  which  formed  as  an  arctic  cold-air  outbreak  and  cold  front  moved  from  the  MIZ 
out  over  the  warmer  ocean. 

Observations  in  the  MIZ  are  limited  by  the  region’s  remoteness  and  the  resulting 
logistic  difficulties  that  this  imposes.  Various  multidisciplinary  experiments  have  made 
observations  of  interactions  between  the  ice,  ocean  and  atmosphere.  The  Arctic  Ice 
Dynamics  Joint  Experiment  (AIDJEX)  was  a  five  year  project  in  the  1970s  to  determine 
ice  pack  motion  and  its  causes.  Using  AIDJEX  data,  Thorndike  and  Colony  (19S2) 
studied  the  causes  of  average  ice  motion.  They  found  the  long-term  average  wind  and 
the  upper-level  circulation  of  the  ocean  equally  contribute  to  the  long-term  (several 
months)  average  ice  motion  for  distances  greater  than  400  km  from  a  land  boundary. 
The  ice  motion  fluctuations  (after  removing  the  mean  motion)  are  more  closely  related 
to  the  gcostrophic  wind  than  to  the  ocean  current.  The  geostrophic  wind  accounted  for 
70%  of  the  variance  of  the  ice  velocity.  Within  400  km  of  land,  the  effect  of  the  coastline 
on  internal  ice  stresses  also  contributes  to  the  ice  motion.  From  the  Norwegian  Remote 
Sensing  Experiment  in  1979,  Johannessen  et  al.  (1983)  found  a  correlation  between  wind 
direction  and  ice  edge  convergence/divergence.  With  on-ice  flow  the  MIZ  is  more 
compact  than  with  office  flow. 

The  most  extensive  documented  ABL  observations  in  the  MIZ  have  come  from  the 
Marginal  Ice  Zone  Expciiment  (MIZEX)  in  the  Bering  and  Greenland  Seas  during  19S3 
and  19S4.  Johannessen  (19S7)  describes  the  summer  MIZEX  East  program  in  Fram 
Strait  and  the  Greenland  Sea.  Fairall  and  Markson  (1987)  show  the  variation  of  surface 
fluxes  for  five  aircraft  flights  during  July  19S3  over  the  MIZ  in  the  Greenland  Sea. 
Conditions  for  the  fi\c  flights  include  winds  from  all  four  directions  relative  to  the  ice 
edge  at  wind  speeds  7  m  s*1  and  below.  They  noted  that  the  wind  direction  and  its  effect 
on  atmospheric  stability,  can  be  an  important  factor  in  describing  the  surface  stress. 
On-ice  flow  or  flow  parallel  to  the  ice  edge  from  the  south  at  less  than  5  m  s*1 ,  caused 
more  stable  conditions  in  the  MIZ.  Under  these  conditions,  the  effects  of  hydrostatic 
stability  on  surface  stress  were  the  same  order  of  magnitude  as  the  effect  due  to  varying 
surface  roughness.  Kellner  et  al.  (1987)  document  turbulence  statistics  and  spectra  un¬ 
der  conditions  of  office  flow  in  the  Fram  Strait.  This  study  includes  aircraft  and 
shipboard  observations  made  during  the  summer  of  1984.  The  mixed  layer  depth  in¬ 
creased  towards  the  ocean  as  convective  mixing  from  the  warmer  ocean  became  impor¬ 
tant.  This  coincided  with  a  shift  in  the  peak  frequency  of  the  vertical  wind  component's 
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variance  spectra  from  1  Hz  (100  m)  over  the  ice  to  A  Hz  (250  m)  over  the  ocean.  Sec¬ 
ondary  low  frequency  peaks  were  evident  over  the  ocean.  These  indicated  secondary 
circulations. 

Early  numerical  models  of  the  MIZ  simulated  off-ice  flow.  Overland  et  al.  (1983) 
used  a  one-layer  slab  model  to  show  the  evolution  of  the  ABL.  The  wind  velocity  in  the 
ABL  and  the  height  of  the  base  of  the  inversion  proved  sensitive  to  the  surface 
roughness  decrease  and  heat  flux  increase  from  ice  to  ocean.  By  50  km  oceanward  of 
the  ice  edge,  the  horizontal  temperature  gradient  had  caused  an  8%  increase  in  wind 
speed.  Roughness  changes  were  responsible  for  a  9%  wind  speed  increase  across  the 
MIZ  (approximately  35  km).  The  base  of  the  inversion  slowly  increased  towards  the 
ocean.  This  study  agreed  with  the  observation  that  an  off-ice  wind  field  contributes  to 
ice  edge  divergence.  Reynolds  (1984)  enhanced  the  slab  model  b>  including  cloud  and 
radiation  effects.  The  addition  of  cloud  and  radiation  effects  for  office  flow,  slightly 
increased  the  slope  of  the  inversion  base  and  slightly  decreased  the  wind  speed.  Such  a 
slab  model  adequately  simulates  the  shallow  boundary  layer  of  off-ice  flow,  but  for  flow 
from  other  directions,  a  more  sophisticated  model  is  needed. 

Bennett  and  Ilunkins  (1986)  used  a  multilevel,  two-dimensional,  higher-order  closure 
model  to  simulate  on-ice  flow.  The  model  included  stratus  cloud  effects  and  sea  ice 
thermodynamics.  They  used  observed  data  to  set  the  model's  initial  conditions  of  wind 
and  potential  temperature  profiles  at  the  ice  edge.  The  processes  of  turbulence,  lateral 
and  vertical  advection  and  radiation  were  simulated.  The  model  showed  good  agreement 
with  the  observed  vertical  soundings  at  various  distances  from  the  ice  edge.  The  base 
of  the  inversion  increased  from  500  m  at  the  ice  edge  to  1000  m  at  150  km  iceward  of  the 
ice  edge.  The  study  confirmed  the  importance  of  including  cloud  top  longwave  cooling 
when  modeling  a  cloud  topped  ABL  under  conditions  of  on-ice  flow.  The  model  results 
suggested  that  the  boundary  laser  cooling  away  from  the  ice  edge  was  mainly  due  to 
cloud  top  longwave  cooling,  not  downward  heat  flux. 

Kantha  and  Mellor  (19S9)  used  a  two-dimensional,  higher-order  closure  model  to 
simulate  off-ice,  on-ice,  and  along-icc  flow.  For  off-ice  flow,  the  model  showed  good 
agreement  with  the  Overland  et  al.  (19S3)  results.  With  on-ice  flow  where  there  was 
strong  heat  flux  from  the  ABL  to  the  ice,  an  inversion  tens  of  meters  in  depth  and  a 
strong  jet  resulted.  Tor  along-icc  flow  with  a  slight  component  either  on  or  off-ice,  a 
strong  thermal  ice  edge  front  with  intense  turbulence  resulted.  The  number  of  detailed 
descriptions  of  the  ABL  over  the  MIZ  from  observations  are  few.  More  observations 
are  needed  tc  validate  the  model  results. 
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B.  CEAREX 

The  Eastern  Arctic  Aircraft  Expedition,  which  took  place  in  March  I9S9,  was  a 
component  of  the  multinational,  interdisciplinary  Coordinated  Eastern  Arctic  Exper¬ 
iment  (CEAREX).  This  expedition  provided  the  opportunity  to  investigate  the  structure 
of  the  ABL  in  the  vicinity  of  the  MIZ  in  the  Barents  Sea  using  the  NOAA  P-3  research 
aircraft. 

Under  conditions  of  flow  nearly  parallel  to  the  ice  edge,  one  flight  was  made  on  24 
March  19S9  to  make  measurements  of  the  ABL,  the  background  atmospheric  structure, 
and  turbulent  fluxes  using  a  gust-probe.  This  paper  discusses  the  observations  of  that 
flight.  This  is  the  first  set  of  observations  to  provide  a  comprehensive  description  of 
ABL  structure  over  the  MIZ  for  flow  nearly  parallel  to  the  ice  edge. 
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II.  THE  DATA 


A.  FLIGHT  DESCRIPTION  AND  CONDITIONS 

The  U.S.  Navy  Fleet  Numerical  Oceanography  Center's  surface  pressure  analysis 
of  1200  Coordinated  Universal  Time  (UTC),  24  March  19S9  (Figure  1)  depicted  two  low 
pressure  centers,  one  west  of  Norway  at  60°  N  and  0°  E  and  one  at  approximately  73° 
N  and  50°  E.  There  was  a  short  inverted  ridge  where  the  isobars  and  wind  trajectories 
curve  south  between  these  two  low  pressure  centers.  In  the  vicinity  of  the  short  inverted 
ridge,  which  included  the  measurement  area  centered  at  approximately  75°  N  and  23° 
E,  the  wind  direction  varied.  The  flow  in  the  measurement  area  was  northeasterly.  To 
the  cast  of  the  measurement  area  the  flow  was  more  northerly,  while  to  the  west  it  was 
easterly. 

During  the  aircraft's  approach  to  the  measurement  area,  stratocumulus  clouds  were 
observed  over  the  open  ocean  with  some  precipitation  beneath  the  clouds.  This  cloud 
layer  thinned  towards  the  ice  edge  and  there  were  no  clouds  over  the  ice. 

The  criterion  for  the  flight  was  that  the  wind  Row  would  be  nearly  parallel  to  the  ice 
edge,  with  a  small  off-ice  component  so  boundary  layer  operations  would  not  be  limited 
by  fog.  With  the  criterion  met,  measurements  were  made  in  the  MIZ  southeast  of 
Spitzbergcn  on  24  March  I9S9  between  74  to  76°  N  and  19  to  27°  E.  The  MIZ  itself 
was  fairly  compact.  The  ice  pack  (with  greater  than  S0%  ocean  cover)  followed  a  10  to 
20  km  zone  of  pancake  ice.  The  ice  edge  was  oriented  northeast-southwest  as  indicated 
by  the  dashed  line  in  Figure  2. 

The  aircraft  speed  during  measurements  was  approximately  100  m  sM  .  The  flight 
pattern  (Figure  3)  included: 

•  Four  200  km  legs  perpendicular  to  the  ice  edge  at  900,  S00, 700  mb  and  45  m  which 
were  almost  exactly  centered  over  the  MIZ. 

•  Two  legs  approximately  50  km  in  length  over  the  ice  at  33  and  150  m,  parallel  to 
the  ice  edge. 


Two  legs  approximately  50  km  in  length  over  the  ice  at  75  and  150  m,  perpendicular 
to  the  ice  edge. 


Figure  1.  Surface,  pressure  analysis 
Source:  Lackmann  et  al.,  (1989) 

•  Two  legs  approximately  50  km  in  length  over  the  ocean  at  45  and  267  m,  parallel 
to  the  ice  edge. 

•  Two  legs  approximately  50  km  in  length  over  the  ocean  at  135  and  267  m,  per¬ 
pendicular  to  the  ice  edge. 

IT  INSTRUMENTATION 

The  ABL  mean  and  turbulent  variables  were  measured  from  the  NOAA  P-3.  Ob¬ 
taining  a  wind  measurement  with  respect  to  the  earth  requires  determining  both  the 
aircraft's  velocity  with  respect  to  the  earth  and  the  velocity  of  the  air  with  respect  to  the 
aircraft.  An  inertial  navigation  system  (INS)  provided  the  aircraft  velocity  and  position 
with  respect  to  the  earth.  The  mean  velocity  of  the  air  with  respect  to  the  aircraft  was 
obtained  using  a  pitot-static  tube.  A  forward  port  on  the  tube  gives  a  total  pressure 
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Figure  2.  NO  A  A- 10  IR  imagery  of  1038  UTC  24  March  1989 

measurement  and  side  ports  give  static  pressure.  The  diiTcrcncc  between  these  two 
pressures  is  a  measure  of  dynamic  pressure.  Combining  this  dynamic  pressure  mcas- 
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urcmcnt  with  ambient  temperature  and  static  pressure  measurements  gives  the  mean 
wind  speed.  The  wind  direction  was  obtained  from  vanes  attached  to  the  tube.  The 
turbulent  wind  measurements  were  obtained  using  a  gust  probe  whose  pcrfonnance  is 
described  by  Grcenluit  and  Gilmer  (1985). 

A  Rosemount  total  temperature  probe  was  used  to  obtain  the  mean  temperature. 
Total  temperature  refers  to  the  ambient  temperature  plus  dynamic  healing  (the  effect 
of  air  flow  on  the  sensor).  A  .025  nun  platinum  wire  sensing  element  experiences 
changes  in  specific  resistance  which  are  proportional  to  the  change  in  temperature.  This 
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same  principal  is  used  by  the  bead  thermistor  which  gave  the  temperature  measurement 
as  part  of  the  turbulence  package.  The  sea  surface  temperature  was  remotely  sensed  by 
a  Barnes  PRT-5  radiometer.  It  compared  the  radiance  of  the  sea  surface  in  the  9.5  to 
11.5  /tm  wavelength  band  with  that  of  an  internal  reference  cavity. 

A  dewpoint  hygrometer  provided  an  absolute  measure  of  humidity  (mass  per  unit 
volume).  A  mirror  is  maintained  at  the  temperature  where  the  moisture  in  the  air  just 
begins  to  condense  on  it.  The  condensation  is  monitored  optically.  This  sensor  was 
slow  to  respond  to  a  change  in  the  measured  variable  compared  to  the  other  sensors. 
It  took  tens  of  seconds  to  reach  equilibrium.  Another  part  of  the  gust  probe  system,  a 
microwave  rcfractomctcr,  measured  the  turbulent  moisture  fluctuations.  This  sensor 
directly  measures  fluctuations  in  the  radio  refractive  index.  Combining  this  measure¬ 
ment  with  those  of  temperature  and  static  pressure  gives  the  water  vapor  density. 

There  was  also  a  lidar  system  on  board  the  NOAA  P-3  which  was  operated  by 
University  of  Washington  personnel.  The  system  transmitted  two  laser  light  beams  ev¬ 
ery  second,  of  wavelengths  .532  and  1.046  jam.  The  transmitted  light  is  strongly  scat¬ 
tered  by  cloud  droplets.  A  telescope  collects  the  returned  light  and  focuses  it  on  a 
detector.  From  there  it  is  amplified,  digitized  every  20  ns  and  recorded.  This  resulted 
in  a  resolution  of  7.5  m.  Cloud  base  and  cloud  top  height  measurements  were  obtained 
using  the  lidar.  Table  I  lists  the  variables,  their  sampling  rates  and  the  sensors  used. 
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Table  1.  INSTRUMENTATION  CARRIED  BY  THE  NOAA  P-3  IN  CEAREX. 


Mean  Variables 

Sampling  Rate  (s"‘) 

Instrument  Source 

u,  V 

1 

Pitot-static  system  com¬ 
bined  with  INS 

T 

1 

Rosemount  Total  Tem¬ 
perature  Probe 

102CH2AF 

Dewpoint 

1 

General  Eastern  120DEP 
Dewpoint  Hygrometer 

Altitude 

1 

Radar  altimeter  (Stewart- 
Warncr  APN159) 

Pressure 

1 

Rosemount  542K2 

Sea  surface  temperature 

1 

Barnes  PRT-5 

Radiometer 

Turbulent  Variables 

Sampling  Rate  (s"‘) 

Instrument  Source 

u,  v,  w 

40 

Gust  probe  combined 
with  INS 

T 

40 

Bead  thermistor 

q 

40 

Microwave  refractometer 

C.  DATA  PROCESSING 

Most  of  the  mean  variable  sensors  were  checked  by  comparing  their  measuiements 
against  routine  observations  taken  at  Bodo,  Norway,  where  the  CEAREX  flights  origi¬ 
nated.  This  comparison  revealed: 

•  The  dewpoint  hygrometer  observations  were  5.6°  too  cold. 

•  Temperature  data  were  within  1°  K  and  pressure  data  within  1  mb. 

The  radiometer  used  to  obtain  the  surface  temperature  was  calibrated  at  NOAA's  Office 
of  Aircraft  Operations  (OAO)  in  Miami,  prior  to  and  after  the  CEAREX  flights.  These 
checks  showed  a  drift  of  2°  K  at  -35°  C  and  1°  K  near  freezing. 
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Inertial  navigation  systems  (INS)  generally  experience  slowly  varying  drifts  with 
time.  It  was  planned  to  use  the  Loran-C  navigational  aid  signal  to  compensate  for  this 
following  the  flight.  However,  tne  Loran-C  signal  was  too  weak  to  use.  So,  the  possi¬ 
bility  exists  that  INS  drift  could  have  introduced  variations  as  large  as  3  m  s*1  in  the 
wind  field  (Shaw,  19S8). 

Greenhut  and  Gilmer  (1985)  described  the  calibration  and  performance  of  the  gust 
probe  system.  The  thermistor  is  calibrated  in  a  wand  tunnel  to  correct  for  dynamic 
heating.  This  correction  is  approximately  4°  C  at  an  air  speed  of  100  m  s_1.  The  char¬ 
acteristics  of  the  microwave  refractometer  are  also  well  known.  The  gust  probe  system 
performance  during  a  set  of  aircraft  maneuvers  (called  the  Lenschow  maneuvers)  is  well 
established.  During  the  CEAREX  flight,  these  Lenschow'  maneuvers  were  executed. 
Later  evaluation  of  the  gust  probe  system  data  during  these  maneuvers  identified  any 
abnormalities  in  the  performance  of  the  system. 

The  gust  probe  system  was  mounted  on  a  boom  attached  to  the  nose  of  the  aircraft. 
The  boom  length  was  chosen  so  its  natural  frequency  is  about  14  Hz.  Frequencies 
greater  than  11.5  Hz  are  filtered  out  of  all  the  gust  probe  system  signals  on-board  the 
aircraft,  prior  to  recording  the  data.  This  avoids  aliasing  the  higher  frequency  natural 
resonances  of  the  boom  and  other  sensors  into  the  recorded  data. 

Two  separate  data  streams  of  mean  and  tuibulent  (gust  probe)  observations  were 
recorded,  with  UTC  being  the  only  common  variable.  The  mean  data  were  initially 
processed  at  NOAA's  Office  of  Aircraft  Operations  (OAO)  in  Miami.  All  the  above 
noted  errors  were  corrected  and  the  data  were  screened  for  unrealistic  values. 

The  turbulence  data  were  initially  processed  at  NOAA's  Environmental  Research 
Laboratories  (ERL)  in  Boulder,  Colorado.  The  turbulent  fluxes  were  determined  using 
the  eddy  correlation  method  at  the  Naval  Postgraduate  School.  It  was  desired  to  show 
the  variation  of  the  boundary  layer  fluxes  compared  to  the  surface  variability  and  also, 
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to  use  Fast  Fourier  Transforms  (FFT)  which  use  less  computer  time  than  other  methods. 
To  accomplish  this,  the  data  (time  series)  were  analyzed  in  102.4  s  segments  corre¬ 
sponding  to  about  10  km  of  sampling  distance.  This  102.4  s  segment  length  met  the  FFT 
requirement  of  2’  data  points  per  segment.  The  sampling  rate  was  40  s'1,  yielding  4096 
(  212 )  points  per  Fourier  transform. 

A 

As  noted  earlier,  the  recorded  gust  probe  system  data  included  frequencies  up  to  1 1.5 
Hz.  The  first  step  was  to  pass  the  time  series  through  a  high-pass  symmetric  compound 
running  mean  subtraction  filter  (Holland,  1968).  This  eliminated  the  unwanted  low  fre¬ 
quency  (mean)  information  from  the  time  series,  leaving  the  high  frequency  (perturba¬ 
tion  or  turbulent)  pail.  The  filter  parameters  were  set  such  that  more  than  90%  of  signal 
variance  was  removed  for  wavelengths  longer  than  10  km  (Shaw  et  al.,  1990).  This  is 
equivalent  to  the  filter  eliminating  90%  of  the  signal  information  for  frequencies  less 
than  .01  Hz.  This  particular  filter  was  chosen  for  its  efficiency  anu  east,  of  implementa¬ 
tion.  Another  benefit  of  using  this  particular  filter  vas  that  its  transfer  function  is 
well-behaved  and  does  not  seriously  distort  the  other  frequencies  m  the  time  series.  The 
filter  passed  a  small  amount  of  signal  variance  for  frequencies  of  .01  to  .03  Hz,  and  all 
signal  variance  for  frequencies  greater  than  .03  Hz.  It  was  desired  to  retain  signal  in¬ 
formation  for  frequencies  of  .01  Hz  and  greater.  The  information  that  was  eliminated 
for  frequencies  .01  to  .03  Hz  was  reinserted  during  a  later  step. 

Next,  Fast  Fourier  Transforms  (FFT)  were  performed  ftn  •.-r.uii  102.4  s  segment  and 
the  spectra  and  cospectra  were  calculated.  The  one-sided  spectra  and  cospectra, 
GC'(j)  and  Gcb(j)  where  a  and  b  are  are  any  two  variables  and  f  is  frequency,  describe  the 
variance  and  covariance  distributed  ever  frequency.  The  variance  and  covariance  are 
defined: 
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The  covariances  represent  turbulent  fluxes  (Sorbjat  1989).  For  example;  the  vertical 
temperature  flux  (w'T')  is  a  covariance  between  the  . "  ufent  or  perturbation  quantities 
^(denoted  by  the  prime)  of  vertical  motion  (w')  and  temperature  ^T'). 

The  third  step  was  to  dividethe  spectra  and  cospectra  by  the  tiansfer  function  of  the 
compound  running  mean  subtraction  filter,  which  was  applied  in  the  time  domain.  The 
transfer  function  describes  the  relation  between  the  signal  input  to  the  filter  and  the  fil¬ 
ter's  output.  It  is  obtained  by  performing  a  Fourier  Transform  on -the  function  de¬ 
scribing  the  filter.  The  covariance  transfer  function  is  a  ratio  between  the  covariance 
after  filtering  and  the  covariance  before  filtering.  Tnc  cospectra  computed  in  step  two 
is  a  ratio  between  co.ariance  after  filtering  and  frequency.  Thus  dividing  the  cospectra 
from  step  two  by  the  transfei  function  gives  the  true  cospectra  (cospectra  before  filter¬ 
ing).  This  relation  is  also  valid  for  the  spectra  if  variance  is  substituted  for  covariance. 
This  step  reinserts  the  .01  to  .03  IIz  information  that  was  removed  by  the  filtering  in  step 
one. 

The  last  step  was  to  integrate  the  true  spectra  and  cospectra  with  respect  to  fre¬ 
quency.  This  gave  the  desired  average  variances  and  covariances  (turbulent  fluxes)  for 
each  102.4  s  segment. 
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III.  MEAN  STRUCTURE 


A.  VERTICAL  STRUCTURE 

Vertical  profiles  were  obtained  from  10  s  averages  during  altitude  changes  between 
the  straight  and  level  measurementdegs.  These  altitude  changes  took  place  at  appr  oxi- 
mately  100  km  from  the  ice  edge  over  the  ice  ai  d  ocean.  The  composite  virtual  potential 
temperature  profile  (Figure  4)  showed  that  the  ABL  was  well-mixed  at  approximately 
252'  K  to  about  150  m  over  the  ice  and  at  approximately  263°  K  to  ab*  at  750  m  ever 
the  ocean.  Above  the  ABL,  the  profiles  over  both  the  ice  and  ocean  display  nearly 
identical  thermal  structure  (+  2*  K).  The  composite  specific  humidity  (q)  profile 
(Figure  4)  had  this  same  general  structure  except  that  over  the  ocean,  the  humidity 
profile  was  not  so  well-mi.  ,ed  in  the  boundary  layer. 

Four  altitude  changes  were  made  over  the  ice.  The  resulting  profiles  individually 
revealed: 


•  Three  of  the  four  virtual  potential  temperature  profiles  were  well-mixed  to  about 
100  m  (Figure  5).  There  was  a  slightly  stable  layer  between  100  and  200  m,  ana  u 
strong  capping  inversion  above  200  m.  The  fourth  profile  which  was  the  farthest 
from  the  ice  edge  differed  in  that  the  well-mixed  layer  was  approximately  70  m  deep 
and  inversion  layer  was  .5°  K  warmer  than  the  other  three  profiles.  This  suggests 
the  inversion  was  sloping  downward  sharply  toward  the  ice. 


•  The  dt  'point  hygrometer  had  a  tendency  to  overshoot  a  step  change.  Two  of  the 
four  available  specific  humidity  (q)  profiles  (Figure  6)  were  obtained  from  altitude 
changes  that  started  within  the  ABL,  These  showed  a  constant  q  value,  charac¬ 
teristic  of  a  well-mixed  layer,  to  a  height  of  150  m.  The  other  two  profiles,  one  of 
which  was  the  profile  farthest  from  the  ice  edge,  were  obtained  from  altitude 
changes  that  started  above  the  ABL.  These  showed  q  variation  within  the  ABL, 
making  identification  of  the  well-mixed  layer  top  less  clear.  However,  for  the  pro¬ 
file  farthest  from  the  ice  edge,  q  began  to  decrease  with  height  at  approximately  70 
m.  This  agrees  with  the  virtual  potential  tempeiature  finding  that  the  ABL  depth 
for  the  profile  farthest  from  the  ice  edge  is  about  one-half  that  found  in  the  other 
profiles. 
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•  The  downstream  profiles  of  virtual  potential  temperature  and  specific  humidity  re¬ 
flected  moistening  of  the  boundary  layer  and  wanning  due  to  surface  heat  flux. 

•  The  wind  speed  profiles  were  seriously  contaminated  by  turns  and  rapid  altitude 
changes.  Using  only  10  s  data  segments  where  the  aircraft  heading  changed  less 
then  .5°s_I  and  the  altitude  changed  less  3  ms-1  ,  the  profiles  (Figure  7)  showed 
that  the  ABL  wind  adjustment  to  gcostrophic  value  occured  almost  entirely  in  the 
inversion.  Between  150  and  250  m,  the  wind  increased  by  8  m  s-1 . 

B.  CROSS-SECTIONS 

Shaw  e:  al.,  (1990)  described  the  structure  of  the  ABL  during  the  flight.  The  data 
from  the  200  km  flight  legs  provided  cross-sections  perpendicular  to  the  ice  edge;  The 
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VIRTURL  POTENTIAL  TEMPERRTURE  OVER  THE  ICE  (K) 


Figure  5. 


Proves  of  virtual  potential  temperature  over  ice 


radiometric  surface  temperature  cross-section  (Figure  S)  had  a  step-like  structure.  Over 
the  ocean,  the  surface  temperature  was  slightly  above  0°  C  at  distances  greater  than  40 
km  from  the  ice  edge  and  just  below  0°  C  within  40  km  of  the  ice  edge.  The  ice  edge 
was  clearly  identified  by  a  marked  decline  in  surface  temperature,  with  leads  indicated 
by  warm  temperature  spikes.  Twenty  kilometers  iceward  of  the  ice  edge,  the  radiometric 
surface  temperature  had  dropped  to  -15°  C  and  by  60  km  to  near  -25°  C. 

Virtual  potential  temperature  traces  (Figure  S)  of  the  900,  800,  and  700  mb  legs  re¬ 
flected  near  horizontal  temperature  contours  which  indicate  there  was  virtually  no 
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SPECIFIC  HUMIDITY  OVER  THE  ICE  (G/KG) 

Figure  6.  Profiles  of  specific  humidity  over  ice 


change  of  temperature  perpendicular  to  the  ice  edge  at  any  altitude  above  the  ABL. 
Therefore,  interactions  with  the  auiface  are  responsible  for  all  cross-ice-edge  thermal 
variability  within  the  ABL.  Ocean-ice  temperature  differences  cause  baroclinity  within 
the  ABL.  Consistent  with  the  small  off-ice  component  to  the  flow*,  the  baroclinity  was 
strongest  in  the  region  40  km  oceanward  of  the  ice  edge. 

Using  the  lidar  pointing  downward  on  the  700  mb  leg  and  upward  on  the  45  m  leg, 
it  was  possible  to  determine  the  top  and  bottom  of  the  stratocumulus  layer.  The 
stratocumulus  laver  (shaded  region  on  Figure  8)  w-as  300  to  400  m  thick  at  100  km 


m 

to 


Figure  7.  Profiles  of  wind  speed  over  ice 

oceanward  of  the  ice  edge.  The  cloud  base  was  at  approximately  225  m.  In  the  vicinity 
of  showers,  the  cloud  base  was  lower.  The  stratocumulus  layer  gradually  thinned  to  no 
clouds  over  the  ice.  The  Jidar-dcrivcd  cloud  layer  base  was  compared  to  the  calculated 
lifting  condensation  level  (LCL).  The  dashed  line  on  Figure  S  identifies  the  LCL.  The 
results  were  generally  consistent,  though  the  LCL  was  approximately  25  m  lower.  The 
cloud  top  was  used  as  an  estimate  of  the  boundary  layer  height  (h).  This  method  could 
not  be  used  to  determine  h  over  the  ice,  since  there  were  no  clouds.  The  potential  tem¬ 
perature  profiles  were  well-mixed  in  the  ABL  and  virtually  identical  above  the  ABL  over 
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the  ice  and  ocean.  The  boundary  layer  height  (h)  was  estimated  by  assuming  that  the 
45  m  leg  potential  temperature  data  represented  the  ABL  well-mixed  value.  Then  the 
height  where  this  same  potential  temperature  value  was  found  on  the  free-atmosphere 
profile  was  considered  the  ABL  depth.  This  method  assumes  a  barotropic  atmosphere 
above  the  inversion  in  the  cross-ice  edge  direction  and  therefore  does  not  allow  for  a 
thermal  wind  component  parallel  to  the  ice  edge.  This  assumption  is  not  quite  correct, 
since  the  influence  of  subsidence  from  the  large  scale  inverted  ridge  would  act  to  increase 
the  temperature  gradient  in  the  sloping  inversion  and  make  a  thermal  wind  component 
parallel  to  the  ice  edge  likely.  However,  the  parameterization  appeared  to  give  a  rea¬ 
sonable  estimate  of  h  over  the  ice.  The  lidar  cloud  top  height  observations  were  in  close 
agreement  with  the  parameterized  h  values  near  the  ice  edge. 

The  specific  humidity  traces  (Figure  9)  showed  that  within  the  ABL,  the  air  over  the 
ocean  was  twice  as  moist  (1.2  g  kg*1)  as  the  air  over  the  ice  (.6  g  kg*1).  The  atmosphere 
above  the  ABL  was  quite  dry  with  specific  humidity  values  of  .4  g  kg*1  just  above  the 
ABL  and  .2  gkg*1  at  approximately  1400  m.  In  agreement  with  the  thermal  structure, 
the  strongest  horizontal  moisture  gradient  was  in  the  region  40  km  oceanward  of  the  ice 
edge. 

Windspecd  (Figure  10)  varied  by  only  I  m  s-1  along  the  200  km  flight  legs  above  the 
ABL.  Within  the  boundary  layer  over  the  ocean,  the  windspeed  was  within  1  m  s*1  of 
the  windspeed  above  the  ABL.  Over  the  ice,  however,  the  windspeed  near  the  surface 
varied  from  11  ms*1  at  the  ice  edge  to  about  6  m s*1  at  100  kilometers  from  the  ice  edge. 
The  strongest  horizontal  windspecd  variation  occured  in  the  region  more  than  40  km 
iceward  of  the  ice  edge. 

The  wind  direction  cross-section  (Figure  1 1)  showed  a  variation  of  approximately 
10°  at  most,  on  all  horizontal  flight  legs.  The  near  surface  wind  was  from  approximately 
25°.  Between  the  surface  and  900  mb  (about  800  m),  the  wind  veered  no  more  than  10 
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Figure  9.  Cross  section  of  specific  humidity:  Flight  path  was  perpendicular  to  the  ice  edge.  The  shaded  area  is 
the  stratocumulus  deck  as  detected  bylidar.  The  dashed  line  over  the  water  is  the  lifting  condensation 
level,  determined  from  45  meter  data.  The  dotted  line  over  the  ice  is  the  estimated  boundary  layer 
height,  also  determined  from  45  meter  data. 


WINOSPEEO  (m/s) 


Figure  10.  Cross  section  of  wind  speed,  perpendicular  to  the  ice  edge. 
Source:  Shaw  et  al.f  (1990) 


WIND  DIRECTION  (deg) 


Figure  1 1.  Cross  section  of  wind  direction,  perpendicular  to  the  ice  edge. 
Source:  Shaw  et  al.,  (1990) 


to  12°  .  Between  altitudes  of  1  and  3  km,  the  wind  backed  20°  .  This  implies  cold 
advection.  However,  some  of  the  variation  may  be  due  to  unaccounted  for  drifts  in  the 
INS  system.  Therefore,  the  wind  direction  changes,  particularily  above  the  ABL  where 
the  changes  are  small,  should  be  regarded  with  some  caution. 

C.  HORIZONTAL  GRADIENTS 

The  flight  pattern  with  its  parallel  and  perpendicular  ice  edge  legs  allowed  horizontal 
gradients  to  be  determined  at  two  levels  over  both  the  ice  and  ocean.  Figure  12  shows 
the  spatial  series  for  u,  v  and  virtual  temperature  and  their  linear  best  fit  from  the  50  km 
segment  at  the  end  of  the  45  m  flight  leg  over  the  ocean,  u  and  v  are  in  ice  edge  coor¬ 
dinates  with  positive  u  towards  the  ocean.  A  linear  best  fit  was  obtained  for  each  vari¬ 
able.  This  provided  a  measure  of  how*  much  a  variable  changed  (its  delta-  value)  along 
the  flight  leg.  The  delta  values  were  used  to  determine  the  components  of  the  horizontal 
gradient  using  a  least  squares  method  for  heights  of  40  and  150  m  over  the  ice  and  45 
and  267  m  over  the  ocean.  Results  are  listed  in  Table  2.  Using  temperature  T  as  an 
example,  values  of  AT,  A.v,  and  Ay  derived  from  each  leg  were  used  in: 

at  dT  a  .  ST  A 
A7=*  -7—  A.v  +  -7—  Av 
ox  oy 

This  is  of  the  form  z  =  ax  +  by.  The  least  squares  method  minimizes  £(z  —  z,)2  with 

i 

respect  to  a  and  b. 


-fa  ^J-z‘  ~  +  M2  =  0 

1 


and 


24 


-Jb  ^[z/  -  (a*/  +  byt)f  =  o 

t 

Expanding  and  rearranging  gives: 


and 


Y 4 Yjpi  -  Jj't'Zs'iyi 


The  gradient  calculation  for  pressure  was  slightly  different.  Figure  13  shows  the 
pressure  series  reflected  the  variation  of  the  aircraft's  altitude.  Using  the  radar  altimeter 
height  (z)  values  and  the  hydrostatic  approximation,  the  pressure  was  adjusted  to  a 
constant  flight  altitude  (z,f/)  then  used  in  the  gradient  calculation. 


8z 


=  ~Pg 
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Figure  13.  Pressure  spatial  series:  from  the  45  meter  flight  leg  over  ocean. 


p  =  pRTv 


where  p  is  pressure,  p  is  density,  g  is  gravity,  R  is  the  specific  gas  constant,  and  Tv  is 
virtual  temperature  which  is  assumed  constant  over  small  altitude  changes.  Substituting 
this  into  the  hydrostatic  approximation  and  integrating  gives: 


where  p„f  is  the  pressure  adjusted  to  a  constant  flight  altitude. 

Figure  14  is  fl-c  specific  humidity  series  from  the  dewpointer  and  the  microwave 
refractometer  (10  s  averages).  The  mean  has  been  subtracted  from  the  rcfractometer 
data.  This  is  why  the  mean  values  differ  between  the  two  series.  The  dewpointer's 
tendency  to  overshoot  a  step  change  affected  the  delta  values.  For  this  reason,  the 
microwave  refractometer  data  were  used  to  calculate  the  horizontal  specific  humidity 
gradients.  Table  2  summarizes  the  horizontal  gradients  over  the  ice  and  ocean.  The 
gradients  are  with  respect  to  the  ice  edge  coordinate  system  where  positive  x  is  towards 
the  ocean  and  positive  y  is  along  the  ice  edge  towards  the  northeast.  The  geostrophic 
wind  was  calculated  using  the  pressure  gradient.  The  other  gradients  were  used  in  de¬ 
termining  the  physical  balances  described  in  the  next  chapter. 
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Specific  humidity  spatial  series:  from  the  d5  meter  flight  leg  over  ocean. 


Table  2.  GRADIENTS  DERIVED  USING  THE  LEAST  SQUARES  METHOD 


Over  ice 

Over  ocean 

40  m 

150  m 

45  m 

267  m 

dp 

(mb  /  km) 

-.00905 

-.02287 

-.02215 

-.01498 

.00364 

.00036 

-.00194 

-.00262 

(K  /  km) 

.02010 

.01547 

.02555 

.03319 

—  (K  /  km) 
oy 

-.00847 

-.00522 

-.02460 

-.02947 

d(} 

(g  /  kg  km) 

.00309 

.00261 

.00089 

.00322 

-.00457 

.00211 

-.00097 

-.00616 

4r~  (m  /  s  km) 
cx 

.04815 

-.0123 

.00604 

.02118 

4f~  (m  /  s  km) 
cy 

.00213 

.01334 

.02745 

.03648 

(m  ( s  km) 
cx  ' 

-.00143 

-.00442 

-.00159 

.00526 

—  (m  /  s  km) 
cy 

.020SS 

.00745 

.01097 

.01119 

D.  GEOSTROPHIC  WIND 

The  gcostrophic  wind  (!')  was  calculated  using  the  gradients  from  the  pressure 
sensor  data. 


=  __L££_ 
Ug  fp  oy 


1  op 
fp  ox 
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Figure  15  shows  the  gcostrophic  wind  over  the  ice  and  ocean.  Vt  at  40,  45,  150  and  267 
m  was  calculated  using  the  pressure  gradients.  Since  900  mb  was  at  all  times  above  the 
ABL,  the  actual  wind  at  that  level  is  considered  gcostrophic.  Over  the  ocean,  the 
gcostrophic  wind  decreased  in  magnitude  between  45  and  267  m.  This  is  consistent  with 
the  ABL  thermal  gradients.  The  thermal  wind  (!',)  defined  by  the  gcostrophic  wind 
change  with  height,  blows  with  low  temperature  to  the  left.  Over  the  ice,  the  gcostiophic 
wind  changes  in  magnitude  and  direction  between  40  and  150  m.  This  is  not  consistent 
with  the  ABL  thermal  gradients  since  V,  defined  by  the  change  of  \\  with  height  docs 
not  blow'  with  low  temperature  to  the  left.  It  is  the  invasion  temperature  giadicnt  noted 
in  Figure  4  that  accounts  for  the  1'  change.  Above  the  ABL,  I'  over  the  ice  and  ocean 
has  comparable  magnitude  and  direction.  Thus,  over  tire  ice,  \  \  is  relatively  constant 
with  height  within  the  ABL  but  changes  significantly  in  the  inversion,  while  over  the 
ocean,  J'  changes  significantly  in  the  ABL.  This  relation  has  implications  concerning 
the  momentum  balance  described  in  the  next  chapter. 
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IV.  TURBULENCE  MEASUREMENTS 


This  section  describes  the  heat,  moisture  and  momentum  balances  in  the  ABL  ap¬ 
plied  at  the  corners  of  the  flight  pattern  over  the  ice  and  ocean.  An  ice  edge  coordinate 
system  is  used  with  positive  x  towards  the  ocean  and  positive  y  along  the  ice  edge  to¬ 
wards  the  northeast. 

A.  HEAT  BALANCE 

The  time  rate  of  change  of  potential  temperature,  neglecting  radiation,  is: 

60  _  60  60  cVO' 

ct  cx  cy  cz 

where  0,  u  and  v  arc  the  mean  potential  temperature  and  velocity,  and  ;•/  and  O'  arc  the 
turbulent  vertical  velocity  and  potential  temperature.  The  coordinate  system  is  such  that 
x  is  perpendicular  to  the  ice  edge  and  positive  toward  the  ocean.  The  advcction  terms 
represent  cross-ice  ( u  and  along-ice  advcction  (v  The  time  rate  of  change  of 

potential  temperature  was  determined  from  two  flight  legs,  along  the  same  course  but 
at  different  altitudes,  which  were  separated  by  approximately  I  hr.  This  is  appropriate 
since  both  flight  legs  were  in  the  ABL,  and  the  potential  temperature  profile  was  well 
mixed.  The  velocity  components,  u  and  v,  were  determined  from  20  km  segments  at  the 
flight  pattern  corners.  The  average  vertical  turbulent  heat  flux  (•••/(?')  was  determined 
(sec  Chapter  II,  sections  B  and  C)  for  the  AO  km  segments  closest  to  the  flight  pattern 
corners.  Tigurc  16  shows  the  v/0'  profiles  and  linear  best  fit  which  was  forced  to  zero 
at  .Sh  where  h  is  the  height  of  the  ABL  (120  m  over  ice,  600  m  over  ocean).  The  best 
fit  was  forced  to  zero  at  .Sii  since  this  is  the  expected  height  where  the  heat  flux  changes 
sign  due  to  entrainment  of  higher  temperature  air  from  above  the  inv  ersion  base. 


* 
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Figure  16.  Vertical  turbulent  heat  flux 

Table  3  shows  the  heat  balance  over  the  ice  and  ocean.  Over  the  ice,  there  was  no 
detectable  change  in  potential  temperature  with  time,  while  over  the  ocean,  there  was 
net  cooling.  Advection  mostly  balanced  the  vertical  heat  flux.  Along-ice-edge  advection 
was  approximately  twice  as  large  as  cross-ice  advection.  The  magnitude  of  the  terms 
was  larger  over  the  ocean  than  ice.  Compared  to  the  other  terms,  the  residual  was  small. 
This  strongly  suggests  that  neglected  effects  did  not  play  a  large  role  in  the  balance. 


* 
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Table  3.  HEAT  BALANCE 


89 

di 

(K  /  hr) 

89 

—  M  ~ 

OX 

(K  /  hr) 

8w'9' 

8z 

(K/hr) 

residual 
(K  /  hr) 

Ice 

0.00 

-0.12 

-0.24 

0.46 

0.10 

Water 

-0.46 

-0.48 

-1.15 

0.91 

0.26 

Source:  Shaw  et  al.,  (1990) 


B.  MOISTURE  BALANCE 

The  time  rate  of  change  of  moisture  is: 

8q  dq  5q  dw'q' 

81  cx  5y  8z 

where  q  is  the  mean  specific  humidity,  and  q'  is  the  fluctuating  part.  The  time  rate  of 
change,  u,  v  and  the  average  vertical  moisture  flux  (w'q')  were  obtained  by  the  same 
method  used  with  the  heat  balance.  Figure  17  shows  the  w'q'  profiles  and  linear  best 
fit  which  was  forced  to  zero  at  h  (150  m  over  ice,  750  m  over  ocean).  The  best  fit  was 
forced  to  zero  at  h  since  upward  moisture  flux  is  expected  to  decrease  with  height  up  to 
the  inversion  base. 

Table  4  shows  the  moisture  balance  over  ice  and  ocean.  The  time  rate  of  change 
reflects  net  drying  of  the  ABL.  Over  the  ocean,  advection  partially  balances  the  vertical 
moisture  flux.  The  magnitude  of  the  residual  indicates  this  description  of  the  moisture 
balance  was  not  as  successful  as  with  the  heat  balance.  The  Figure  17  best  fit  over  the 
ocean  does  not  reflect  the  observed  increase  of  the  moisture  flux  with  height  over  the 
lower  ABL.  This  increase  with  height  indicates  drying  near  the  surface  and  moistening 
of  higher  levels.  This  suggests  there  is  a  moisture  source  in  the  ABL  other  than  the  flux 
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Figure  17.  Vertical  turbulent  moisture  flux 

from  the  ocean  surface.  A  similar  increase  of  w'q'  with  height  over  the  ocean  was  noted 
in  a  study  by  Nicholls  and  Readings  (1979).  Their  observations  showed  w'q'  increased 
with  height  over  the  lower  ABL  and  decreased  with  height  in  the  upper  ABL.  The  mean 
specific  humidity  decreased  with  height  throughout  the  mixed  layer.  They  concluded  the 
increasing  w'q'  profile  was  acting  to  change  the  decreasing  mean  profile  to  a  mor  well 
mixed  (constant)  profile.  The  situation  during  the  CEAREX  flight  differed  slightly  in 
that  the  mean  specific  humidity  profile  didn't  decrease  with  height  throughout  the  mixed 
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layer.  However,  Figure  4  shows  it  was  not  so  well  mixed  at  the  lower  levels.  No  ob¬ 
servations  were  available  to  determine  if  w'q'  decreased  with  height  in  the  upper  ABL. 

The  moisture  flux  increase  with  height  implies  a  downward  transfer  of  water  vapor 
in  the  boundary  layer.  For  this  to  occur,  a  moisture  source  above  the  surface  is  required. 
Donelan  and  Mivake  (1973)  suggested  evaporation  of  water  droplets  released  from  the 
surface  by  whitecapping  as  a  possible  moisture  source  that  could  explain  the  flux  in¬ 
crease  with  height. 


Table  4.  MOISTURE  BALANCE 


Gq_ 

Ol 

(g  /  kg  hr) 

8q 

—  u  — — 

ox 

(g  /  kg  hr) 

1 

Gw'q' 

Gz 

(g  /  kg  hr) 

residual 
(g  /  kg  hr)  ■ 

Ice 

-.001709 

-.00612 

-.001235 

-.0078S 

.01353 

Water 

-.0564 

-.0137 

-.14799 

.29124 

-.18595 

C.  MOMENTUM  BALANCE 

The  momentum  balance,  neglecting  accelerations  is: 


Gu'vS 

Gz 


Gv'v/ 

Gz 


=f[u-  ug) 


(1) 

(2) 


The  above  equations  sa>  that  the  difference  between  the  pressure  gradient  force  (fvs)  and 
the  coriolis  force  (/V)  is  balanced  by  the  momentum  flux  divergence  (stress).  Therefore, 
the  difference  between  the  observed  and  geostrophic  wind  is  related  to  the  stress  profiles. 
If  the  difference  is  constant,  the  stress  changes  linearly  with  height.  If  the  difference  isn't 
constant,  curvature  in  the  stress  profile  results.  The  observed  wind  averages  for  the 
flight  levels  within  the  ABL  arc  listed  in  Table  5. 
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Table  5.  OBSERVED  WIND  AVERAGES 


height 

u  (m/s) 

v  (m/s) 

Ice  75m 

3.05 

-6.78 

Ice  40m 

1.47 

-6.75 

Ocean  267m 

3.65 

-12.22 

Ocean  135m 

3.75 

-11.67 

Ocean  45m 

4.8 

-11.95 

1.  Over  the  Ice 

Table  5  indicates  that  u  varies  and  v  is  constant  with  height  over  the  ice.  Recall 
that  over  the  ice,  Vg  did  not  vary  with  height  within  the  ABL  since  the  boundary  layer 
was  too  shallow  for  the  thermal  wind  to  have  a  significant  effect.  So,  when  applying  the 
momentum  balance  in  the  ABL,  Vt  is  assumed  to  be  constant  with  height.  Since  vs  and 
v  are  considered  constant,  their  difference  is  constant  and  equation  (1)  requires  the 
profile  of  u'vS  to  be  linear.  Figure  18  shows  the  mV  profile  and  linear  best  fit.  The 
coordinate  system  is  oriented  such  that  the  mV  profile  is  perpendicular  to  the  ice  edge. 

The  u  wind  component  varies  in  the  boundary  layer.  Assuming  that  u  varies 
linearly  with  height  allows  u  to  be  described  by: 

>  Am 

U  =  U°  +  ~KTZ  (3) 

where  u„  is  u  at  the  surface  and  z  is  the  height  above  the  surface.  Integrating  equation 
(2)  over  height  gives: 


rV(z)-vV(0)  =  f  *M-u)dz 


(4) 
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Figure  18.  Vertical  turbulent  momentum  flux  over  the  ice 


Substituting  equation  (3)  in  equation  (4): 


(5) 


Integrating  and  rearranging  yields: 


v'\v'{z)  =  vV(0)  +Jz  ug  -  u{ 0) -y—j-z 


(6) 
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Equation  (6)  describes  the  best  fit  curve.  Using  the  v'w'  averages  from  the  various  alti¬ 
tudes,  a  least  squares  method  was  used  with  this  curve  to  obtain  v'w'(O).  Table  5  values 
were  used  for  ,  and  equation  (3)  was  solved  for  u(0).  Figure  18  shows  the  v'w' 
profile  and  best  fit  defined  by  equation  (6).  The  coordinate  system  is  oriented  such  that 
the  v'w'  profile  is  parallel  to  the  ice  edge.  The  curve  closely  matches  the  observed  v'w' 
averages  of  the  flight  legs  perpendicular  to  the  ice  edge  while  there's  some  scatter  for  the 
flight  legs  parallel  to  the  ice  edge.  The  closeness  of  the  best  fits  to  the  observed  values 
confirms  that  the  neglected  inertia  terms  did  not  have  a  dominant  role. 

2.  Over  the  Ocean 

The  large  heat  flux  and  surface  stress  produced  a  strong  mixed  layer  with  little 
change  of  observed  winds  with  height  (Table  5).  Recall  that  over  the  ocean,  the  mag¬ 
nitude  of  Vt  varied  with  height  in  the  ABL  (Figure  15)  due  to  the  temperature  gradient. 
Table  6  lists  the  geostrophic  wind  components  over  the  ocean. 


Table  6.  GEOSTROPHIC  WIND  AVERAGES 


height 

ut  (m/s) 

v,  (m/s) 

Ocean  267m 

1.44 

<N 

CO 

i 

Ocean  45m 

1.04 

-11.90 

Table  6  indicates  that  ut  is  relatively  constant  and  vt  varies  with  height.  Since  ut  and  u 
are  considered  constant,  their  difference  is  constant  and  equation  (2)  requires  the  profile 
of  v'w'  to  be  linear.  Figure  19  shows  the  v'w'  profile  and  linear  best  fit. 

Assuming  that  v?  varies  linearly  with  height  allows  v,  to  be  described  by: 


Av„ 


y  si*  4-  ■ 

*  S°  +  Az 


(7) 


where  vi0  is  vt  at  the  surface.  Integrating  equation  (I)  over  height  gives: 
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Figure  19.  Vertical  turbulent  momentum  flux  .over  the  ocean 


(8) 


Substituting  equation  (7)  in  equation  (8): 


uV(z)  -  uV(0)  = 


(9) 
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u v(z)  =  uV(O)  +fz  V  -  vg(0)  -  y  2 

Equation  (10)  describes  the  expected  functional  form  of  the  variation  of  wV  with  height. 

Av, 

Table  6  values  were  used  for  »  an^  equation  (7)  was  solved  for  v4(0).  Using  the 
u'w'  averages  from  the  various  altitudes,  a  least  squares  method  was  used  with  this  curve 
to  obtain  u'w'{ 0).  Then  u'w'(0)  was  used  in  equation  (10)  to  obtain  the  best  fit. 
Figure  19  shows  the  u'w'  profile  and  best  fit.  The  negative  u'w'  of  the  45  m  flight  leg 
remains  unexplained.  There  were  no  indications  of  equipment  malfunction. 
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V.  SUMMARY 


Data  collected  with  the  NOAA  P-3  research  aircraft  was  used  to  examine  the  mean 
structure  and  physical  balances  of  the  ABL  in  the  MIZ.  The  wind  flow  was  nearly  par¬ 
allel  to  the  ice  edge  with  the  900  mb  geostrophic  wind  approximately  15°  off-ice.  The 
atmospheric  conditions  above  the  ABL,  were  virtually  identical  at  100  km  from  the  ice 
edge  over  both  the  ice  and  ocean.  Therefore,  the  structure  of  the  ABL  was  due  to 
interaction  with  the  surface.  The  well  mixed  boundary  layer  extended  at  least  100  km 
into  the  ice.  The  horizontal  gradients  of  temperature  and  moisture  were  maximum  at 
approximately  40  km  oceanward  of  the  ice  edge  but  no  boundary  layer  front  was  evident. 

The  heat  balance  showed  that  the  vertical  heat  flux  was  mostly  balanced  by 
advection.  Over  the  ocean,  the  magnitude  of  both  the  heat  flux  and  advection  was  quite 
large.  The  residuals  were  small  compared  to  the  other  terms.  This  strongly  suggests  that 
neglected  effects  did  not  play  a  major  role. 

The  moisture  balance  gave  an  unusual  result.  Over  the  ocean,  the  profile  of  vertical 
moisture  flux  increased  with  height  over  the  lower  ABL.  This  profile  could  be  viewed 
as  acting  to  make  the  mean  specific  humidity  profile  more  constant  or  it  could  be  de¬ 
scribing  the  presence  of  a  counter-gradient  flux. 

The  relationship  between  the  observed  and  geostrophic  winds  in  the  momentum 
balance  explained  the  curvature  in  the  stress  profiles.  The  boundary  layer  was  suffi¬ 
ciently  deep  over«the  ocean,  to  allow  the  horizontal  temperature  gradient  to  affect  the 
geostrophic  wind.  Over  the  ice  the  ABL  was  shallow,  thus  minimizing  this  effect. 

These  observations  provide  evidence  of  the  interactions  occuring  in  the  MIZ  under 
conditions  of  flow  nearly  parallel  to  the  ice  edge.  The  sloping  inversion  and  horizontal 
temperature  gradient  had  a  significant  effect  on  the  geostrophic  wind  in  the  ABL.  The 
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geostrophic  forcing  in  turn  affected  the  surface  stress.  The  magnitude  and  variability 
of  these  interactions  are  important  for  understanding  ice,  ocean  and  atmospheric  dy¬ 
namics  in  the  MIZ. 
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